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' Abstract 
O ■ 

, We derive an elementary formula for Janossy densities for determinantal point processes with 

_ ' a finite rank projection-type kernel. In particular, for (3 = 2 polynomial ensembles of random 
matrices we show that the Janossy densities on an interval / C M can be expressed in terms of 

' the Christoffel-Darboux kernel for the orthogonal polynomials on the complement of /. 

in 

^ ; 1 Introduction 

\ We consider an ensemble of n particles on a measure space (X, //) with the joint distribution density 

' (with respect to the product measure /i®") given by the formula 

(N 



(N 
O 



^ \ where X" = X x • • • x X (n times). Ensembles of this form were introduced in j3Uj . In the special 

' cases when X = (j)i = ipi = x^~^, and X = G C| |z| = 1}, 0j = = z^~^, such ensembles 

^ were extensively studied in Random Matrix Theory much earlier under the general name unitary 

ensembles, see |21| for details. An example of the form which is different from random matrix 
ensembles was considered in |22j . 

Let us assume that and {V'j}i=i,...,n with respect to the pairing 

= / <l)ix)Tp{x)fi{dx). 
Jx 
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p{xi, ... ,Xn) = COnstn ■ dct ((/^j (xfc) )j,fc=l,...,n dct (V'j (xA;))j,fc=l,...,n • (1) 

Here (j)k{x), ipk{x), k = 1, . . . ,n, are some functions on X and constn is the normalization constant 



const„ ^ 



/ det{<j)j{xk))j^k=i,...,ndet{'ipj{xk))j,k=i,...,n TT fJ-idxj] 



g : = n! det ( / (l)i{x)ijj{x)fi{dx) ) (2) 



X / j,j=l,...,n 
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In other words, suppose that we can find functions £,k{x), rik{x), k = 1, . . . ,n such that 

^fc e Span(0j, j = 1, . . . ,n), % € Span(V'j, J = 1, • • • , n), {^k,Vrn) = ^km- 

The famihes {^j} and {rjj} are called biorthogonal bases in Span((/)j, j = 1, . . . ,n) and Span(V'j, j = 
1, . . . ,n) considered as subspaces in ^). Then the distribution can be rewritten as (P, pUj l 

. . . = ^ det{K{xi,Xj))ij=i^,„^n, (3) 

with 

n 

K{x,y) = Y,^ji^H{y)- (4) 

One of the particularly nice properties of the ensemble © is that one can explicitly calculate 
the correlation functions 



Pk{xi, . . . ,Xk) := TTT / p{xi,...,Xk,Xk+l,...,Xn)lJ.{dXk+l)--- IJ.{dXn) 

[n-ky. Jx^^-k 

which still have a determinantal form with the same kernel K{x,y): 

Pk{xi, ... ,xk) = det{K{xi,Xj))ij=i^,„k. (5) 

If fj, is supported by a discrete set of points then the probabilistic meaning of the /c-point corre- 
lation function is that of the probability to find a particle at each of k sites xi, X2, ■ ■ ■ Xk- In other 
words, 

Pk{xi,X2, ■ ■ ■ Xk)p{xi) ■ ■ ■ fj,{xk) = Pr {there is a particle at each of the points Xi, i = 1, . . . , k}. 
Analogously, if X C M and fi is absolutely continuous with respect to the Lebesgue measure then 
Pk{xi,X2, . . . Xk)P'{dxi) ■ ■ ■ p{dxk) = Pr {there is a particle in each interval {xi,Xi + dxi)}. 

In general, random point processes with the /c-point correlation functions of the determinantal 
form © are called determinantal or fermion (see e.g. _271). 

So-called Janossy densities Jk,i{xi^ ■ ■ ■ , Xk), k = 0,1,2, ... , describe the distribution of the par- 
ticles in a subset I X. If X C M and p is absolutely continuous with respect to the Lebesgue 
measure then 

Jk,i{xi, ■ ■ ■ Xk)p{dxi) ■ ■ ■ p{dxk) = Pr {there are exactly k particles in I, 

one in each of the k infinitesimal intervals {xi,Xi + dxi)}. 

If p is discrete then 

Jkj{xi, . . . Xk)p{xi) ■ ■ ■ p{xk) = Pr {there are exactly k particles in I, 

one at each of the k points Xi\. 

See m] for details. 
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For determinantal point processes Janossy densities also have a determinantal form (see |12j . 
p. 140, or j4j, section 2): 

Jkj{xi,.. . = const{I) ■ <lei{Li{xi,Xj))i^j=i^,„k, (6) 

where 

Li = Kj{Id-Ki)-\ (7) 

Here the kernel of Kj is the restriction of the kernel K{x,y) to /: Ki(x,y) = y)x/(y), 
where xi{ ' ) is the characteristic function of /, and const{I) is the Fredholm determinant 

const{I) = det(/d - Kj) = det{Id + Lj)'^. 

The main result of this paper is 

Theorem 1. Let ^j, j = l,...,n and rjj, j = l,...,n be biorthonormal bases in Span{(pj, j = 
1, . . . , n} and Span{'ipj, j = 1, . . . ,n} considered as subspaces of \ /, /j,): 



^k & Span{(t>j, j = 1,. . . ,n), 7]k e Spaniipj, j = 1, . . . ,n), / ^k{x)rim{x)fi{dx) = 6km- 

Jx\i 

Then the kernel of Lj = Kj(Id — Ki)~^ is equal to 

n 

Lj{x,y) =Y,U^)Vjiy)- (8) 

i=i 

The above result readily applies to the so-called f3 = 2 polynomial ensembles. Such ensembles 
arise, in particular, in random matrix theory |2^, |10| . directed percolation and tiling models |18| 
ElUni) and representation theory j^l, [HI- The definition is as follows. 

Assume that X is a subset of M and take 

p{xi,...,Xn) =COnstn- {Xj-Xkf. (9) 

(Recall that this formula gives the joint distribution density with respect to ^{dxi) ■ ■ ■ fj,{dxn)-) 

This is a special case of with 4'j{x) = ipjix) = x^~^, j = 1, . . . ,n. Then we have = rjj = 
Pj-i, where {pj{x)} are normalized orthogonal polynomials on {X,fi{dx)), and deg{pj) = j. The 
kernel K{x,y) is the nth Christoffel-Darboux kernel 

K.(x,y) = X;p,(x)p,(,) = b^ Pn[x)p^-M-Vn{y)Vn-.{x) 

where kj is the coefficient of x^ in Pj{x). It should be noted that the kernel K depends on n, but in 
what follows we will usually omit the subscript n unless this may lead to a confusion. 

Clearly, in the case of the polynomial ensemble @, Theorem states that the kernel of Lj = 
Kj{l — Kj)~^ is the nth Christoffel-Darboux kernel computed for the measure ^ restricted to X \ /. 
That is, 

kn-lPn{x)Pn-l{y) -Pn{y)Pn-l{x) 



j=o ^-y 
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where 

Pj{x) = kjX^ + {lower order terms}, / pk{x)pm{x) ii{dx) = 5km- 

Jx\i 

One of the particulary nice properties of the Janossy densities is that for any interval (or, more 
genertally, a measurable set) / and non- negative integer k one has 



Pr( there are exactly k particles in /) = — / J'k i(xi, 

k\ J Ik ' 



,Xk)fJ.{dxi) ■ ■■ fj,{dxk) 



The Janossy densities can be particularly useful in calculating the distribution of the left-most 
(right-most) particles when the particle space X is a subset of the real line. Indeed, let us denote by 
-^1 !^ ^2 < . . . < A„ the locations of the particles in the increasing order. Then it is easy to see that 

Pr(Afc e {s,s + ds)) = i ^ I JkA~oo,s)ixi, • • • , Xk-i, s)fx{dxi) ■ ■ ■ fi{dxk-i) ] /i(ds) 

\(^ - 1)' J(-oo,s)fe-i J 

= Pr(Ai > s) ^ I det(L(„oo,s)(a;i,x-,))ij=i,...fc^((ixi) • • •/Li(dxfc-i)^((is) (10) 

Pr(Ai >s) = (det(Id + L(_^,,)))"^ = det(/d - i^(_oo,s)), 

(where in l\W\i we put Xk = s). 

This observation and the Theorem above allow us to compute explicitly the distribution functions 
of the left-most particles in the hard-edge scaling limit of random matrix models when the parameter 
(charge at the edge) is equal to zero. We refer to Section 4 below for the details. 

The result of Theorem 1 was initially discovered in the case of polynomial ensembles using the 
techniques of Riemann-Hilbert problems. Later on, it was realized that Theorem 1 has a simpler 
linear algebraic proof. However, since the main idea of the "Riemann-Hilbert" computation is very 
useful in deriving Painleve equations for the distribution of the left- or right-most particles in de- 
terminantal point processes, see (Zj, (21, |S1) we decided to include the argument into this paper; it 
can be found in Section 2. Section 3 contains the simpler proof. Concluding remarks are given in 
Section 5. 

To conclude the Introduction, let us note that Theorem 1 has a counterpart for the so-called 
pfaffian ensembles. (The (3=1 and 4 (or "orthogonal" and "symplectic" ) random matrix ensembles 
are the most known examples of the pfaffian ensembles.) See the companion paper j2H] for details. 



2 Riemann-Hilbert Problem 

In this section we will briefly describe two applications of the Riemann-Hilbert problem (to comput- 
ing orthogonal polynomials and to inverting integrable integral operators) and use them to derive 
Theorem 1 in the case of polynomials ensembles. Since we use the Riemann-Hilbert problem (RHP, 
for short) mainly for instructional purposes, we avoid the discussion of any technical issues involved. 

Let S be an oriented contour in C. We agree that when we go along the contour in the direction 
of orientation, the positive side lies to the left and the negative side lies to the right. Let vhe a, map 
from S to GL(Z, C), where Z = 1, 2, . . . . We say that an / x / matrix function m = m{z) is a solution 
of the RHP (S,w) if (0, [ini) 

(i) m{z) is analytic in C\S, (11) 
(m) m+(z) = m_(z)f (z), z G S. (12) 
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Here ■m^(z), m-{z) stand for the limiting values of m{z) as z approaches S from the positive 
(negative) side. If, in addition, m{z) Id as z ^ oo then we say that m{z) solves the normalized 
RHP {T,,v). The matrix v is usually called the jump matrix for the RHP. 

First we describe the connection of RHP to orthogonal polynomials, see ^Ij, |15j . Let djji[x) = 
Lo{x)dx be an absolutely continuous measure on the real line, such that the non-negative density co 
decays at infinity sufficiently fast (in particular, all moments exist). Consider an RHP on M oriented 
from left to right with the jump matrix 

= ( ; '^<;' ) . (13) 

Fix a non-negative integer n. We are looking for a solution of the RHP (M, v) satisfying 



i{z) = {ld + 0{z~^)) 





m[z) = [ia + U[z ^)) \ ^ I, z^oo. 



It appears that this RHP has a unique solution given by 

m{z)=( ^^rr^^^lv^V ^^^^ (14) 

V 7n-i7r„_i(z) 7„_i(C(a;7r„_i))(z) J 

where 7r„(z) = z" -|- . . . is the nth monic orthogonal polynomial corresponding to the weight function 
a;(x), 

(CA)u)=J-/m,f 

2vrz Jni- z 

is the Cauchy transform, 7„ = — 27rifc^, and /c„ is the leading coefficient of the nth orthonormal 
polynomial p„, i.e. Pn{z) = knTTniz). Thus, computing the orthogonal polynomials with the weight 
uj{z) is equivalent to solving RHP of the form above. 

Now let us explain the relation of RHP to integrable operators. Let / be a subset of M (typically, 
a disjoint union of finitely many intervals). We recall that an integral operator M in L^(/, dx) with 
the kernel M{x,y) is called integrable (HSl, |211, HH) if 

M(^^^y) = lLiM^Myl (15) 

x-y 

for some / = 2,3, . . . and some functions /j, gi on /. We assume that Yli=i fi{x)9i{^) = so that 
the kernel has no singularity on the diagonal. 

In particular, the formula for the Christoffel-Darboux kernel in the case of polynomial ensembles 
discussed in §1 means that the operators K and Kj can be viewed as integrable operators in L^(R, dx) 
and L'^{I, dx) with / = 2, and we may take 

fi{x) = ^;^^J^Pn{x)^/uJ(x), f2{x) = -kn-iPn~i{x)\/uj{x), (16) 

gi{x) = 27ri/c„„ip„_i(x)-\/u;(x), 52(2;) = Pn{x)y^Lo{x). (17) 



The appearance of y'uAx) has to do with the fact that we consider the Lebesgue measure rather 
than fj,{dx) as our reference measure for the L^-space. 
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It turns out that if the operator Id — M is invertible then the resolvent R = M{Id — M) ^ is 
also an integrable operator and 

i?(x,y) = S^i^lM^, (18) 
x-y 

Fi = ild-M)~^f,, Gi = {Id-M')-^g,, i = l,2,...,/, (19) 

see pi] and also |llj for a very nice exposition. Furthermore, the functions Fi and Gi can be 
obtained through solving a RHP as follows. Let be an / x / matrix valued function on I given by 

v' = Id-2mfg\ f={f,,...JiY, g={g^,...,gi)\ (20) 

One can prove ([23]) CB) that the normalized RHP {I,v') has a unique solution m'{z), and 

F = (Fi,...,FO* = (m')±/, (21) 
G = (Gi,...,Q)* = (m')i*5. (22) 

The following observation is crucial. 

Lemma 1. Let m be the solution Ji^l ) of the RHP (R, v) and let m' be the solution of the normalized 
RHP (I, v') with v' given by ^EW and (l^. Then M = m'm solves the RHP (M \ /, v) with the 

asymptotics diag{z^ , z~'^) as z ^ oo, and hence 

Miz) = ( ^ , , ^ (7^V^!!))(^) , y . R \ /, (23) 

V ln-iT^n-i{z) 7,„„i(C(u;xR\/7r„_i))(z) J 

where ~ signifies that the corresponding polynomials are orthogonal on M \ / with respect to the same 
weight function uo . 

The proof of this lemma is based on Lemma 4.3 of [7] (see also Lemma 2.4 of "7 for a discrete 
analog). The analog of Lemma 1 for weight functions with discrete support was one of the basic 
tools used in 0. 

Proof. A straightforward calculation shows that on I we have v' = mj^v~^m^ = m-v~^niZ'^ . Thus 
on / 

MZ = ml {m!_)~ m'^nij^ = mZ v'tuj^ = v'^ mZ = Id. 

On the other hand, since m'{z) is holomorphic away from / and tends to Id as z ^ oo, it is clear that 
on M\ /, M(z) satisfies the same jump condition as m{z), and that it also has the same asymptotics 
as m{z) when z — > oo. □ 

Proof of Theorem 1. We apply the formalism described above to the Christoffel-Darboux kernel with 
/j, gi specialized by (fTC]) . (fT7|) above. As before, we will use the notation p^, tt/j for the fcth orthonor- 
mal and monic orthogonal polynomials corresponding to the weight w on M \ J, and we also denote 

gfc = C(w7rfc), qk = C{LJXM.\fik)■ 
\a the calculations below we use the identity det m(2;) = detM(z) = 1. Indeed, Liouville's 
theorem readily implies that if the jump matrix of a RHP has determinant 1 and the determinant 
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of the asymptotics of a solution at infinity is also equal to 1, then the determinant of any solution of 
this RHP (having the corresponding asymptotics at infinity) must equal 1 identically. 
We have 

^1 = {m'f)i = (Mm"^)ii/i + (Mm~^)i2/2 = Miim22/i - Mi2m2i/i - Miimi2/2 + Mi2mii/2 

n—lPnQnPn—1 knqnPnPn— 

= -^^{QnPn-1 - Qn-lPn)PnV(^ = det(m)j3„V^ = ~ PnV^- 

knkn ZTTlkn 2mkn 

Similar calculations yield 

Gi = 2mkn-lPn-iy/^-, F2 = -kn-lPn-lV^, G2 = :^PnV^. 

kn 

Hence, the kernel of L/ = Kj{Id — Kj)~^ equals 

Fi{x)Gi{y) + F2{x)G2{y) kn-l Pn{x)Pn-l{y) - Pn{y)Pn-l{x) 



Li{x,y) = = -~ ^/uJ{x)u;{y). 

x-y kn ^-y 

Recall that the factor ^/uj{x)io{y) is due to the fact that we are working in L'^{I,dx) rather than 
L'^{I,uj{x)dx). The proof of Theorem 1 for polynomial ensembles is complete. 

□ 



3 Linear algebraic proof 

We use the notation of §1. Consider the integral operators Kj and Lj in //) with the kernels 

n 

Kn{x,y) =J2^j{x)rij{y) (24) 

and 

n 

Li{x,y) = Y,U^)m{y)- (25) 

Both operators are finite-dimensional: 

Ran(i^'/) = Ran(L/) = Tii := Span(0j) = Span(^j) = Span(^j), 
K.e£{Ki) = Ker(L7) = ?^^, 7^2 := Span(y;j) = Span(r7j) = Span(rjij), 

where the index j ranges over {1, . . . , n}, and the subspaces are taken inside L^(/, jj). Therefore, in 

order to prove that L/ = Kj{Id — Kj)~^, it is enough to prove this relation for the restrictions of 
Kj and Lj to the n-dimensional space Hi. To this end we compute the matrices of the restrictions 
of the operators Kj, Lj to Hi in the basis {Cj}j=i,...n- 

Let us denote by Gj and Gx\i the n x n matrices with entries 

{Gi)jk= I ij{x)j]k{x)n{dx), {Gx\i)jk= i ij{x)if]k{x)n{dx) j,k = l,...,n. 
Ji Jx\i 
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Since {i]k} are biorthonormal on X, we have Gj + Gx\i = Id. The matrix of Kj on TCi in the 

basis {^j} is given by Gj. To calculate the matrix of the restriction of Lj on TCi we biorthonormalize 
the functions ^j, r]j, j = 1, . . . n, in L'^{X \ I, fi). This gives (cf. Proposition 2.2 in J^) 

n n n 

It immediately follows from ()26|) that the matrix of the restriction of L/ to in the basis {^j} is 
equal to Gj{Id — Gi)~^. The proof is complete. 

4 Hard edge with zero charge 

In a few special cases, the polynomials orthogonal with respect to u; on X \ I can be easily expressed 
in terms of the orthogonal polynomials on X. We consider the Laguerre ensemble of positive definite 
matrices as an example. 

Every positive definite nxn matrix M can be written (in a non-unique way) as M = AA* , where 
A is an n X n matrix with complex entries and A* is the adjoint matrix. The probability measure in 
the Laguerre ensemble (also called Wishart ensemble in statistics) is defined as ((HI): 

P{dM) = Z-^ exp(-Tr(ylyl*)) det{AA*)"dA, (27) 

where dA is the Lebesgue measure on the 2n^-dimensional space of n x n complex matrices, 

is a normalization constant, and a > — 1. The joint probability density of the distribution of the 

eigenvalues of M = AA* is equal to 

n 

p{xi,. . . ,Xn) = constn (xj - Xj)^ a;"e~^^ Xj £ {0, +oo), j = I, . . . ,n. (28) 

l<j<j<n j=l 

The polynomials orthogonal with the weight lu{x) = x"e~^ on M_|_ = (0, -|-cxd) are the classical 
Laguerre polynomials (see e.g. JSl)- In the special case of a = and I = (0, i), the orthogonal 
polynomials on X \ I = M_|_ \ (0, t) = [t, +oo) are obtained from the Laguerre polynomials by the 
simple shift of variable x x — t, i.e. Pj{x) = pj{x — t), j = 0, 1, ... . Thus, by Theorem 1, the 
kernel of lJ'"^^") = K,^"^^"^ (jd - K^^3{n)y^ ^^^^j 

Lj'"^("^(x,y) = -t,y-t), (29) 

where K^^^^"\x,y) is the order n Christoffel-Darboux kernel for Laguerre polynomials with a = 0. 

It is well known, see [2^1, ISHI, that when n becomes large, the smallest eigenvalues in the 
Laguerre ensemble are of order n~^. Moreover, if we rescale all the eigenvalues of the nth Laguerre 
ensemble by then there exists a scaling limit as n ^ oo of all the correlation functions. (In the 
random matrix theory this procedure is usually referred to as "hard edge scaling limit".) The limit 
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correlation functions also have the determinantal form ((S)) with the so-called Bessel kernel: 
hm n-V^'("'"^ (^,...,^) =det(K(")(x„x,)),,=i, 



n n 



»1 oo 
/ Ja+l{2^/Tx)Ja+l{2y/Ty)dT = y2 

Jo 



x-y 



^ kT{a + A; + 1) /!r(a + / + a + k + l + l' 

Here Jv{') is the J-Bessel function, see e.g. jl3j . Note that if a = then the above formula makes 
sense for any rr, y S C. 

Proposition 1. For any s > 0, let K^''^ be the (bounded) integral operator in L^((0, s), dx) defined 
by the restriction of the Bessel kernel K^"'^(x,y) with a = to (0,s) x (0, s). Then the operator 

Ks^^ (l — Ks^^j is bounded and has a kernel which is equal to K^'^^x — s,y — s). 



Proof. The relation p9)) implies 



Since n ^K^"'3^'^\xn ^,yn -'^) tends to -fC^''-' (a;, y) as n —> oo uniformly on compact subsets of C (this 
follows, e.g., from the proof of Theorem 4.5 in taking the scaling limit in (|3()|) yields 

K^^\x,y) = K^^\x-s,y-s)- f K^^\x - s,u - s)K^^\u,y)du, x,ye{0,s). (31) 



Denote by L^^^ the operator in L^((0, s), dx) with the kernel L^s^\x,y) = K^^\x — s,y — s). Since 
the kernel of this operator is the uniform limit of the kernels of the nonnegative operators i/^"'^^'^^ = 
— t,y — t), we have L^^^ > 0, and hence —1 does not belong to the spectrum of L^^\ Thus, 

H31() can be rewritten in the form Kg^'' = L^^^ {^d + L^'^^^ , and this is equivalent to the statement 
of the proposition. □ 

Corollary 1. Let A^"^ < Ag"^ < • • • < An"^ be the ordered eigenvalues of the Laguerre ensemble \2^) 
with Q = 0. Then 

Pr (a^") > ^) = e-^ (32) 
lim Pr (a^i > -) = / det(i^(°)(x„ Xj)),j=i,...,fc dxi • ■ ■ dxk, k>2. (33) 

In particular, 

lim Pr ( A^"^ > = ^ r'^x(/2(x) - lf{x))dx, (34) 

n^oo \ nJ 2 Jq 



where is the I-Bessel function. 
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The formula (|32|) was first observed in jl6j . The hmiting distribution (|34|) of the second smallest 
eigenvalue was computed in ^] and |17j . ^ Further results in this direction, including formulas 
similar to H33|) can be found in j32j . 

Proof. The relation is easy: 

/(«/n,+oo)" n (2^* - n e-^^dxj /(o,+oo)n n ix^ - xj)^ n e-^^-^/"dxj 

Kj j Kj j 

The relation H33() follows from (|lUp applied to the Laguerre ensemble and the uniform convergence of 
kernels mentioned in the proof of the proposition above. Finally, using the L'Hopital rule we obtain 

J\(^){-x,-x)dx = ({J'o{2iV^f - '^^(^'^^^(^'^^ _ Jo(2iV^)J^'(2iV^)) dx. 

The formulas 

Joiz) = -Jiiz), Jq{z) = z'^Ji{z) - Jq{z), Io{z) = Joiiz), Ii{z) = -iJi{iz), 
see |13| . and the change of variable x ^ x^/2 bring the last integral to the form ()34() . □ 

Of course, if Ai < A2 < • • • are the ordered particles of the determinantal point process with the 
correlation functions given by the Bessel kernel with a = then the right-hand sides of (|32|) and 
are equal to Pr (Ai > s) and Pr (A^+i > s), respectively. 

The calculations similar to those above can be done for the Jacobi ensemble corresponding to 
a;(x) = (1 X G (-1, 1), in the special cases / = {t, 1), a = 0; / = f3 = 0. After 

appropriate rescaling one again obtains the limit relations of the form (|32|) , (|33)) , (|34j) . 



5 Concluding remarks 

Take two n-point ensembles with joint probability densities of the form const-pn{xi, . . . , Xn) with the 
same pn, but assume that these two ensembles are supported by different sets — the first one lives on 
{X, n) while the second one lives on {X\I, where / is a subset of X. Of course, the normalization 
constants for these two ensembles will be different. The kth Janossy density Jk,i{xi, ■ ■ ■ x^) of the 
first ensemble is given by the formula 

Jkj{xi, ... Xk) = const' ■ / . . . ,a;n)/i((ixfc+i) • • • /i(dxn), xi, . . . , Xfc G /, 

while the kih. correlation function pk{xi, . . . , Xk) of the second ensemble equals 

Pk{xi, ... ,Xk) = const" ■ / Pn{xi, . . . ,Xn)n{dXk+l) ■ ■ ■ p{dXn), Xi , . . . , Xfc G X \ /. 

J{X\I)'^-'' 

The only difference between the two formulas above is in the constant prefactor and in the domain 
where allowed to vary. Clearly, this suggests that there should be a direct relation 

^Note that the integral l|34|l can be evaluated in terms of Bessel functions, see e.g. (2.30) in |29| . 



10 



between Jk i and pk, and in the case of determinantal ensembles such relation is provided by Theorem 

1.2 

Since the argument above does not depend on the specific form of the density pn, one might 
expect that Theorem 1 should have an analog for the pfaffian ensembles (see e.g. |2H], [201, [Hl| for 
definitions). This is exactly the case, and the corresponding result is presented in the companion 
paper |28] . 
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